Large detuning limit for the multipartite systems interacting with electromagnetic 

fields 
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We study the dynamics of the multipartite systems nonresonantly interacting with electromag- 
netic fields and show that, since the coupling strength g is collectively enhanced by the square root 
of the number of microparticles involved, the more rigorous large detuning condition for neglecting 
the rapidly oscillating terms for the effective Hamiltonian should be A > y/N g, instead of A ^> g 
usually used in the literature even in the case of multipartite systems, with A the detuning and 
and N the number of microparticles. This result is significant since merely the satisfaction of the 
original condition will result in the invalidity of the effective Hamiltonian and the errors of the 
parameters associated with the detuning in the multipartite case. 
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I. INTRODUCTION 

Microparticles or artificial microparticles interacting 
with electromagnetic fields are important coupling mod- 
els to accomplish various quantum information process- 
ing tasks, such as quantum gate [1—3], quantum memory 
[4, 5] and quantum communication [6, 7]. Many efforts 
have been devoted to the investigation of these physical 
processes theoretically [8-10]. 

Recently, James and Jerke presented a useful compact 
formula for deriving an effective Hamiltonian describing 
the time-averaged dynamics of detuned quantum systems 
[8]. Then Camel and James developed a technique for 
finding the dynamical evolution of the time-averaged den- 
sity matrix, and the equation of evolution includes an 
effective Hamiltonian, as well as decoherence terms in 
Lindblad form [9]. Brion et al. derived an unambigu- 
ous effective two-level Hamiltonian of a lambda system 
excited by off-resonant laser beams and quantified the ac- 
curacy of the approximation achieved [10]. In this paper, 
we focus on the large detuning condition for neglecting 
the rapidly oscillating terms for the effective Hamilto- 
nians of the multipartite systems interacting with the 
electromagnetic fields. 

In the large detuning regime, the energy exchanges 
between the field modes and the relevant microparti- 
cle transitions are suppressed and the rapidly oscillating 
terms in the original Hamiltonian could be neglected ap- 
proximately to get a simpler effective Hamiltonian, which 
no longer contains the energy exchange terms between 
the field modes and the relevant microparticlc transitions 
accordingly. The large detuning condition is usually for- 
mulated as A > 3 even in the case of multipartite sys- 
tems, actually, which is sufficient only for few micropar- 
ticles involved. When a large number of microparticles 
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interact with the fields, the coupling strength g could be 
collectively enhanced due to the many-particle interfer- 
ence effects [6, 11], therefore the large detuning condition 
in multipartite situation should scale with the number of 
microparticles involved accordingly. We show that the 
rigorous large detuning condition should be expressed as 
A ^> \^Ng in the multipartite case. In fact, if merely 
A g is satisfied, there may exist notable real energy 
exchanges between the field modes and the relevant mi- 
croparticlc transitions in multipartite case, however the 
effective Hamiltonian contains no energy exchange terms 
between them, meaning its invalidity of describing the 
coupled system. In mathematical language, A ^> g is 
necessary but not sufficient for the validity of the effec- 
tive Hamiltonians of the multipartite systems interacting 
with fields, while A 3> \/~N g is just necessary and suffi- 
cient for it. Although by the full numerical simulations 
for this type of problem, one may observe that the detun- 
ing should be chosen much larger than the value deter- 
mined by the condition A» j, this is a time-consuming 
task [12]. Furthermore, for the theoretical works without 
full numerical simulations, the invalidity of the effective 
Hamiltonians resulting from the original large detuning 
condition may not be realized [13—16], and the incorrect 
detuning value will lead to other wrong parameters asso- 
ciated with it, such as the gating time. 



In what follows, we take the cavity QED system as ex- 
ample to derive the large detuning condition for the effec- 
tive Hamiltonian of many atoms interacting with fields. 
First we calculate the exact solution of the time-evolution 
operator for the two-level atoms nonresonantly interact- 
ing with the cavity mode and numerically simulate the 
state evolution. Then we analyze the cases of three-level 
A atoms off-resonantly interacting with the cavity mode 
and the classical field under the conditions of two-photon 
resonance and nonresonacc. 
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II. TWO-LEVEL CONFIGURATION 

First, let us consider an interacting system consisting 
of an ensemble of atoms with two- level configuration non- 
resonantly coupled to a single-mode cavity field. The 
atoms are assumed for simplicity to have the same cou- 
pling strength g and detuning A. In the interaction pic- 
ture, the coupled system can be described by the Hamil- 
tonian (assuming h = I) 
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where the subscript j represents the jth atom, a (a^) 
denotes the annihilation (creation) operator for the cav- 
ity mode, TV is the number of atoms involved and o~ = 



\9 3 )(e. 



> J j 



\ej)(9j\ with \ej), \gj) (j 



1,2, 



being the excited and ground states of the jth atom. For 
suppressing the cavity decay, we require no real energy 
exchange between the cavity mode and the atomic en- 
semble, and the cavity mode in vacuum state initially. 
To this end the coupled system need to work in a large 
detuning regime usually expressed in the form of A ^> g 
in literatures even in the case of a large number of atoms 
involved [13, 15]. While we will show this is not sufficient 
to prevent the real energy exchanging between the cavity 
and the multi-atomic system. 

In the context of quantum information processing with 
microparticlc ensembles [1, 6, 17, 18], one usually con- 
siders the singly excited situation, i.e., only one atom 
is excited in the ensemble while the other atoms in 
the ground states and the cavity mode in the vacuum 
state, or the cavity mode is in the single-photon state 
while all the atoms in the ground states. In this case 
there exists an invariant subspace spanned by {|+) = 
|-E)|0),|— ) = |G)|1)} for the system evolution, where 

\E) = (l/VN)Y;?= 1 \9i---9 j -ie j 9 j+ i---9N), \G) = 
lffi.92 ' ' ' 9n)i an d |0) (|1)) denotes the cavity mode in the 
vacuum state (in the single-photon state). We may call 
\E) the collective atomic excited state and \G) the col- 
lective atomic ground state [4, 19]. In terms of collective 
atomic operators, the Hamiltonian (1) can be written as 
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where S = \G)(E\ and Sft = \E)(G\ are collective atomic 
lowering and raising operators. We can see that the cou- 
pling constant is enhanced by the square root of the num- 
ber of atoms involved, however the detuning is still A. 
From this Hamiltonian, as well as the collective atomic 
level diagram shown in Fig. 1, if we expect there is no 
real energy exchange between the cavity mode and the 
atomic ensemble, then A ^ VNg is required for the 
large detuning regime, instead of the original condition 
A g. In the following, we give a rigorous derivation 
for this result through solving the exact solution of the 
temporal evolution of the coupled system. 
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FIG. 1: The collective atomic level diagram for the two-level 
configuration, (a) Considering the singly excited situation. 
The single-atom-cavity coupling constant is g and the detun- 
ing between single atomic resonance and cavity mode is A. 
(b) The coupling strength between the collective atomic states 
and the cavity mode is y/N times larger than the single-atom- 
cavity coupling constant, however the detuning between the 
collective atomic resonance and the cavity mode is still A. 
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FIG. 2: Temporal evolution of the state |+) under atom 
and cavity decays with the initial state in |+). The three 
solid curves from top to bottom correspond to the number of 
atoms in the cavity being 1, 5, and 25, respectively. Other 
parameters: A = Wg, k = g/10, and 7 = g/100, where 
k is the cavity decay rate and 7 is the atomic spontaneous 
emission rate. 



After a rotating transformation with respect to 
-AEf=i^f/2 with a] = |e,-)(e,-| - \ 9j ){gj\, the Hamil- 
tonian (1) goes into 



JV 



Hr 



(3) 



In the subspace spanned by {|+) = 175)10),!-) = l^)! 1 )}* 
the Hamiltonian Hr can be written as 



Hr = Hi + H 2 , 



with 
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where / is the 2x2 identity matrix. The time-evolution 
operator of the system can be written as 



U(t) = e 



-itH le -itH 2 



(0) 
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Neglecting the global phase factor resulting from Hi, 
through direct calculation the exact solution of the time- 
evolution operator U (t) is given by 



U(t) 



K Bp 
Bp A+ 
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where 



A~ = cos(iaA<) - ^sin(iaAt), 
A+ = cos(iaAi) + £ sin(ioAi), 
Bp = —i/3 sin (|aAi) , 
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In the large detuning regime the effective Hamiltonian is 
given by [!)] 
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where A = g 2 /A. From Eq. (9), we see that there is no 
real energy exchange between the cavity mode and the 
atoms, thus the off-diagonal elements Bp in the tempo- 
ral evolution operator should approximate to zero all the 
time, and according to Eqs. (8) we obtain A ^> VNg, 
in accordance with the physical picture shown in Fig. f . 
With numerical simulation employing master equation, 
figure 2 shows the time evolution of the state |+) under 
atom and cavity decays assuming the initial state in |+). 
We can see that, when the number of atoms in the cavity 
grows, the system is not only more likely to go out of the 
state |+) but also decays faster, which both indicate ex- 
plicitly that the original large detuning condition become 
inefficient gradually with atom number increasing. 



III. THREE-LEVEL A CONFIGURATION 




FIG. 3: The collective atomic level diagram for the A-level 
configuration, (a) Considering the singly excited situation. 
The single-atom-cavity coupling constant is g with detuning 
A + 8 and the coupling strength between the single atom and 
the classical field is with detuning A. (b) The coupling 
strength between the cavity mode and atoms is enhanced by 
the square root of the number of atoms involved, while the 
detuning is still A + 5. 



Consider the singly excited situation similar as in the 
case of two-level configuration, then the system evolu- 
tion will be restricted in the Hilbert space spanned by 
{\m) = \E)\0), \u 2 ) = \G)\1), K) = |5>|0)}, where 

\S) = {1/y/N) Ylj=i 1 5i • ■■9j-i s i9j+i ■ --9n)- In terms 
of collective atomic operators, the Hamiltonian (10) can 
be written as 
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s — \E) (S\. We can see that 
the effective coupling strength between the atoms and 
the cavity mode is y/~N times larger than the single-atom- 
cavity coupling constant, while the detuning between the 
\G) <-> \E) transition and the cavity mode is still A + 6 
(as shown in Fig. 3). Then the excited state |e) of short 
lifetime could be adiabatically eliminated during the op- 
erations, if the excited state is not occupied initially and 
A ^> y/Ng,£l (A 3> g is not sufficient). In what fol- 
lows, we prove this result in two situations: two-photon 
resonance and nonresonance, and give the condition for 
eliminating the cavity mode in the latter case. 



Now, let us consider an ensemble of atoms with three- 
level A configuration (shown in Fig. 3) simultaneously 
interacting with a single-mode cavity and driven by a 
classical field. The two lower states \g) and \s) can be 
achieved, for example, with hyperfine or Zeeman sub- 
levels of electronic ground states, and |e) is the excited 
state. The \g) f-> |e) transition is nonrcsonantly cou- 
pled to the cavity mode with the atom-cavity coupling 
strength g and detuning A + 6, while the \s) •<->• \e) tran- 
sition is nonresonantly driven by the classical field with 
the Rabi frequency f2 and detuning A. In the interac- 
tion picture, the coupled system can be described by the 
Hamiltonian (assuming K = 1) 
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A. Two-photon resonance situation 

First, let us consider the simple case of two-photon 
resonance, i.e., 6 = 0. After a rotating transformation 
[10] with respect to -A^]^ |e,-) (e,-|, the Hamiltonian 
(10) goes into 



N 



H R = £ [A \ ej 

3=1 



{g\ ej )(g 3 \a + n\ 



Sjl+H.c.)]. 
(12) 

In the subspace spanned by {\ui) = \E) |0), \u 2 ) = 
\G) |1), 1^3) = I if?) |0)}, the Hamiltonian Hr can be writ- 
ten as 
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FIG. 4: (Color online) Temporal evolution of the state = 
\E) |0) with the initial state in \ui) = \G) |1) considering no 
dissipations. With the number of atoms in the cavity growing, 
there exists more real energy exchanges between the cavity 
mode and \g) <B> |e) transition in a certain detuning A. Other 
parameters: Q. — g and A = lOg only enough for few atoms 
involved to prevent the \g) <H» |e) transition. 



The exact solution of the temporal evolution operator up 
to a global phase factor e ~** A / 2 is given as follows: 
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A, = T^j2 elAt ^ 2 + i + i/ v 2 A+, 
A/,-T+W ejAt/2 + TlV^' f] = VNg/n, 

A± = cos {\Aat) ± £ sin (\Aat) , 
B/j = -/3sin (|Aat) , S 7 = -7 sin (|Aai) , 
a = v/4iV. 9 2 + 4f7 2 + A 2 /A, 
/3 = 2ViV.g/ v /4iV. g 2 + 4ft 2 + A 2 , 
7 = 2VL/^/ANg 2 +4,n 2 + A 2 . 
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In the large detuning regime, we can eliminate the ex- 
cited state adiabatically, that means the transitions be- 
tween the lower states and the excited state are forbid- 
den, thus Bp and B 7 should approximate to zero all the 
time and according to Eqs. (15) A ^> y/Ng,Q should be 
satisfied. Fig. 4 shows that, as the number of atoms in- 
volved grows, the excited states not occupied initially are 
populated with higher probability during the state evolu- 
tion in a certain detuning A, which reveals that the large 
detuning condition should be relevant to the number of 
atoms involved. 

It should be pointed out that there may exist real en- 
ergy exchanges between the cavity mode and atoms even 
in the large enough detuning due to the two-photon res- 
onance, however it is between the cavity mode and the 
\g) <H- \s) transition, not between the cavity mode and 
the \g) O |e) transition. In the following, we analyze the 



two-photon nonresonance situation, where if the detun- 
ing 5 is large enough, the real energy exchanges between 
the cavity mode and the \g) <-> \s) transition could also 
be blocked. 



B. Two-photon nonresonance situation 



In this case (5 7^ 0), we perform a rotating transforma- 



tion with respect to Y^,jLi [(^ + A) \gj) 
the Hamiltonian (10) goes into 
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In the subspace spanned by {|ui) = \E) |0), \u 2 ) = 
\G) |1), \u 3 ) = \S) |0)}, the Hamiltonian Hr can be writ- 
ten as 



Hr = H 1 + H 2 , 



with 
Hi 



-N (5 + A) /, H 2 = 
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where / is the 3x3 identity matrix, and Hi only re- 
sults in a global phase factor, thus we could neglect it 
in the derivation of the state evolution. Assuming the 
initial state \ip(t = 0)) = a \u±) + f3 \u 2 ) + 70 K3), the 
derivation of the exact solution of the state evolution 
can be performed by finding the eigenenergies of the 
system and using the initial condition to determine the 
coefficients of the Fourier decompositions of the differ- 
ent amplitudes. Since the exact solution of the state 
evolution is too complicated to write down and ana- 
lyze, we introduce a set of reduced variables (A, A C ,A;) 
such that Ae = S/A,X c e = \^Ng/A,Xie = fi/A, with 
< e < 1, (A, A c , A;) = O(l), to simplify the expres- 
sion of the state evolution but preserving the most im- 
portant components. We do not reproduce the calcula- 
tions but only give the final result of the state evolution 
\tp(t)) = a(t)\ Ul ) + (3{t) \u 2 ) + 7 (t) |ua) , with 
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where {pi,p 2 ,P3} are the solutions of the equation 
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and the coefficients {a.,-, /3j, jj }j=\ t 2 3 are determined by 
the initial condition. The expansions in e of these differ- 
ent parameters are given by 

pi =1 + Ae + O (e 2 ) , p 2 = O (e 2 ) , p 3 = Xe + O (e 2 ) , 
ai =a + (A)A C + 7 A;) e + O (e 2 ) , 
"2 = - /?oA c e + O (e 2 ) , 
Q!3 = - 7oA;e + O (e 2 ) , 
ft =a A c e + O (e 2 ) , 
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Assuming the initial state is \u2) = \G) |1), i.e., (ao = 
0, (3q = 1,70 = 0), then the solution of the state evolution 
\<p(t)) = a(t) \ui) + (i{t) \u 2 ) + "f(t) IM3) can be obtained 
from the coefficients: 



a{t) =0 (e 2 ) e- ltAp3 + 
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In this case if we expect there exist no transitions be- 
tween the lower states and the excited state, and no 
real energy exchange between the cavity mode and the 
\g) <-> \s) transition, then both a(t) and 7 (t) should ap- 
proximate to zero all the time, thus it should be satisfied 
that A > -/Ng, 5 > VNgtt/A. In fact, VNgtt/A is 
the collective Raman Rabi frequency, which should be 
much smaller than the Raman detuning 6 for prevent- 
ing the real energy exchanges between the cavity mode 
and the Raman transition. Similarly, assuming the initial 
state is I7/3) = \S) |0) we can get the condition A ^> fl 
for preventing the transitions between the states |s) and 
|e). Fig. 5 shows that, as the number of atoms involved 
grows, the Raman transitions between the states \g) and 
\s) happen with higher probability during the state evo- 
lution in a certain detuning <5, meaning that there will be 
notable real energy exchanges between the cavity mode 
and atoms in the case of large number of atoms involved. 
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FIG. 5: (Color online) Temporal evolution of the state = 
\E) |0), |ua) = \G) |1), |r*3> = \S) |0)} with the initial state in 
^2) = \G) |1) considering no dissipations. The solid curves 
correspond to the state |«3} = \S) |0), the dashed curves cor- 
respond to the state [112) = \G) |1) and the dotted curves 
correspond to the state |ui) = \E) |0). With the number of 
atoms in the cavity growing, there exists more real energy ex- 
changes between the cavity mode and \g) +-> \s) transition in 
a certain detuning 5. Due to the large detuning A the prob- 
ability of the state |wi) = \E) |0) is very small all the time. 
Other parameters: = lOg, A = lOOp and S = 0.3g. 



IV. DISCUSSION AND CONCLUSION 

The enhanced coupling strength between the photons 
and atoms is due to the many-atom interference effects, 
and actually the many-photon interference effects also 
exist as the atoms coupled to the field mode. In fact, the 
coupling strength between the photons and atoms could 
be enhanced not only by the square root of the num- 
ber of atoms, but also by the square root of the number 
of photons involved, thus the large detuning condition 
is also relevant to the photon number if there are more 
than one photon participating in the atom-photon inter- 
acting. In addition, though all the analyses above are 
based on the cavity QED system with atoms, the con- 
clusion is also applicable to other kinds of microparticles 
and artificial microparticles interacting with fields, such 
as ions, polar molecules and superconducting devices, for 
the cavity model and free space model. 

In conclusion, we have studied the dynamics of the 
multipartite systems off-resonantly interacting with elec- 
tromagnetic fields, focusing on the large detuning condi- 
tion for the effective Hamiltonian in this case. Two-level 
configuration and three-level lambda configuration under 
the conditions of two-photon resonance and nonresonacc 
are analyzed theoretically and numerically in detail, and 
the exact solutions of the temporal evolution of the cou- 
pled systems are derived. Since the coupling strength 
is enhanced due to the many-atom interference effects, 
we suggest a more rigorous large detuning condition for 
the effective Hamiltonian in this case. It is significant to 
apply the rigorous large detuning condition in analyzing 
the multipartite systems interacting with fields, since the 
original large detuning condition may lead to the inva- 
lidity of the effective Hamiltonian and the errors of the 
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parameters associated with the detuning, especially for 
the analyses without full numerical simulations. 
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